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Radiation constraints from cosmic strings 
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We show that it is possible to evolve a network of global strings numerically including the effects of radiative 
backreaction, using the renormalised equations for the Kalb-Ramond action. We calculate radiative corrections 
to the equations of motion and deduce the effect on a network of global strings. We also discuss the implications 
of this work for the cosmological axion density. 



1. INTRODUCTION 

Cosmic strings can play many interesting roles 
in the early universe. It has been shown that a 
network of cosmic strings will evolve towards a 
self-similar scaling solution on the largest scales. 
However, on the smallest scales, where radiative 
effects are important, the evolution is less well 
understood. Here we report on an ongoing re- 
search project[l-3] to understand the true nature 
of these effects. Our primary motivation is to 
calculate accurately the amplitude and spectrum 
of the radiation background from both local and 
global strings. For the purposes of this paper we 
shall concentrate on Goldstone boson radiation 
from global strings, however analogous calcula- 
tions have been made for gravitational radiation 
from local strings[4]. 

2. RADIATIVE BACKREACTION 

The starting point for our discussion of global 
strings will be the Kalb-Ramond action[7] 
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which has been shown to exhibit the same dynam- 
ics as the U(l) Goldstone model[l] using com- 
parisons between numerical field theory simula- 
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tions and analytic calculations of the radiation 
power. The antisymmetric tensor field is 
dual to the massless Goldstone boson associated 
with the string, -ff^a/j is the field strength asso- 
ciated with B^^ and rfcr'"' = e"'' daXt' diX" is the 
area element in terms of the worldsheet coordi- 
nates X''((T, r) = (r, X((T, r)). 

Varying the action with respect to the world- 
sheet coordinates and the antisymmetric tensor 
field yields the string equations of motion and 
field equations. In the temporal transverse string 



gauge, X 







t 



T and X • X' = 0, the equations 



of motion for the worldsheet coordinates are 
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where 

F" = {f°,ef + f°X) = 27rUH''"^Vap, (3) 
and 

X'2 



1-X2 
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The field equation can be simplified using the 
analogue of the Lorentz gauge in electromag- 
netism, 9^5'"' = 0. Using this gauge, the field 
equations are 



□5„/3 = 27r/a y" dadf6^{x-X(a,f))Vap. 



(5) 



This equation can be solved easily using stan- 
dard retarded Green function techniques. In anal- 
ogy to electromagnetism, one can perform a series 
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of manipulations to deduce the Lienard-Wiechart 
potential[9], 
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and its derivative, 
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where = a;^ - X^((T,r), A?^^^ = and 
tr < t. Using the sum and difference of the re- 
tarded and advanced time Green functions, one 
can calculate the self- and radiation-fields respec- 
tively. However, since the Lienard-Wiechart po- 
tential and its derivative both involve an integra- 
tion along the string, which in the case of long 
strings is infinite, we must make some approxima- 
tion in order to calculate the field at some point 
on the string. We shall employ the 'local backre- 
action approximation '[3]. In this approximation 
the integrals in (6) and (7) can be truncated at 
some scale A, which at present is arbitrary. One 
can then expand the expressions as a power se- 
ries in A/R (assuming A < R), where R is the 
average radius of curvature of the string. 

Retaining leading order terms, one can deduce 
expressions for the forcing terms in the equation 
of motion (2) as 
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where /u(A) = /jq + 27r/^ log(A/(5) is the renor- 
malised string tension. The finite radiation back- 
reaction force experienced by the string due to 
£rad ^j^j jO,rad jg ^j^g analogue of the Abraham- 
Lorentz force for the point electron in classical 
electrodynamics. 

Performing numerical simulations of the point 
electron using the Abraham-Lorentz force is haz- 
ardous, since the triple derivative with respect to 
time leads to unphysical exponentially growing 
solutions. However, in the case of strings one can 
make the approximation X = 'X."/e^, which re- 
duces the backreaction force to the analogue of a 
viscosity term[3]. Once one has made this approx- 
imation, methods similar to those used in ref.[6] 
can be used to evolve general string trajectories. 

For general long string perturbations parame- 
terized by wavelength L and relative amplitude 
e = 2'itA/ L, where A is the actual amplitude, the 
radiation power per unit length due to the back- 
reaction force can be seen to be 
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(13) 



for some constant j3 ^ f^. This leads to exponen- 
tial decay of e, since the energy per unit length in 
the perturbation is proportional to . This was 
verified by numerical simulations of string tra- 
jectories. We also compared the decay of string 
trajectories to those seen in numerical field the- 
ory. The rate of decay of amplitude was found 
to agree, if the correlations in the long range 
fields were suppressed at the curvature radius and 
A ~ R/A. We also observed that an initially 
sharp kink visibly rounds under the action of this 
force. 
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where 6 ^ f^^ is the string width. 

The self-force can be clearly seen to be a mul- 
tiple of the equations of motion [8, 9], which facil- 
itates the well-known renormalisation, 
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3. NETWORK EVOLUTION 

Our current understanding of the evolution of 
a cosmic string network is based on a marriage 
between analytic models and sophisticated net- 
work simulations [5]. However, the network sim- 
ulations just evolve the free equations of motion 
for a string in an expanding universe. In order to 
incorporate the radiative effects discussed in the 
preceding section one must modify the equations 
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of motion to include a radiation damping term, 
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where a is the scale factor. 

However, to calculate this radiation damping 
term one must use Green functions in an expand- 
ing background[10]. In the radiation era, the re- 
tarded Green function is 

Dretix, X') = (^(^ _ ^')2^ - ^' ) . (16) 

where = (?y,x) and x'f' = (ry',x'). Applying, 
the 'local backreaction approximation' in this sce- 
nario, one finds that the forcing terms are given 

by, 
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where the fiat suffix denotes the fiat space back- 
reaction force given by (10) and (11), and (g°,g) 
is a correction to the force due to the expanding 
background. 

The forced equations of motion (15) can be 
used to derive equations for the evolution of the 
density of long strings (poo) and loops (ph), under 
the infiuence of the expansion, Hubble damping 
and the radiation backreaction force. If one now 
inserts a term to take into account of loop pro- 
duction, the equations become 



Poo = (1 + (V ))Poo f- 
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(v ) is the average string velocity, do,di,.. are 
constants and c is a measure of the efficiency of 
loop production. 



If one now substitutes po 



pC/t' 



id L 



C'^I'H into (19), then one can deduce a differen- 
tial equation for C,. This equation has an attrac- 
tive fixed point, which corresponds to the scaling 
regime. If di = for i > 0, then one can deduce 
that 
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In the case where di is non-zero. One should ob- 
serve transient effects in the scaling. However, for 
large times these effects will be negligible and the 
attractive fixed point is exactly that for di = Q. 

4. RADIATION CONSTRAINTS 

The evolution of a network of global strings is of 
particular interest, since the best motivated sce- 
nario for the production of global strings is at the 
Peccei-Quinn phase transition. In this case, the 
radiated Goldstone boson is the initially massless 
axion, which acquires mass at the QCD phase 
transition, and may form a substantial part of the 
dark matter in the universe. In the light of the 
discussion of the previous section, we concluded 
that a network of global strings will evolve by 
the production of loops, which subsequently radi- 
ate into Goldstone bosons[12]. This mechanism is 
similar to that generally accepted to be the case 
for local strings. However, since the radiation into 
Goldstone bosons is about three orders of magni- 
tude stronger than that into gravitational radi- 
ation, the loop production size a and the string 
scaling density ( are likely to be different. 

Using the scaling balance arguments of §3, one 
can deduce that the loop number density is given 

by[2] 



n(l,t) 
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where k is the loop radiation backreaction scale. 
From this one can deduce that the spectral den- 
sity of axions emitted by the loops is 
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Assuming axion number conservation from the 
QCD phase transition to the present day one can 
deduce that the contribution to the density of the 
universe from axions produced by the loops is 
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where To is the current temperature of CMBR, 
the Hubble constant is Hq = lOO/i km s~^Mpc~^ , 
A ~ 10^*^'® takes into account the uncertainties 
of the QCD phase transition and O.l^a/K^l.O. 
This contribution is substantially larger than any 
of the other possible contributions from long 
string radiation and domain wall collapse. It is 
also much larger than original estimates which as- 
sumed axion production via coherent oscillations 
of a scalar field. This difference can be under- 
stood intuitively as being due to the renormalised 
string energy per unit length. 

5. DISCUSSION 

We have summarized the current understand- 
ing of string network evolution and its relevance 
to constraints on the axion. The formalism we 
described in §2 offers the prospect of evolving 
a string network numerically, with inclusion of 
radiative effects. We have shown that the ef- 
fect of the radiation backreaction force on the 
largest scales of the string network will be sim- 
ilar to that in the most recent analytic work[ll]. 
However, since the Goldstone boson radiation is 
much stronger than gravitational radiation, real- 
istic network simulations of global strings should 
see the scaling of small scale structure much ear- 
lier than for local strings. 

At present the most accurate estimate of the 
axion density due to string loops (25) is depen- 
dent on the ratio a/n. This ratio expresses the 
key uncertainty in string network evolution and 
should in principle be calculable, to within a fac- 
tor of two, in the near future. 

Assuming A 1, To 2.7K and imposing the 
constraint Sla < 1, one can deduce a constraint 
on /a. If a/K 1 then 



/a<4.4 X 10^°GeV, h = 1.0 , 
whereas if a/n 0.1 then 
/a<2.5 X 10^°GeV, /i = 0.5 , 

/a<8.0 X lO^^GeV, h = l.O. 



(27) 

(28) 
(29) 



At their most extreme, the parameter uncertain- 
ties A/}"2[To/(2.7K)]^ can vary in the range 0.2 
to 25, therefore for O.l^a/K^l.O the constraint 
must lie in the region 



/a<10^ - 10"GeV, 



(30) 



Note that even the weakest string bound is 
stronger than the early estimate due to coherent 
zero-momentum axions, /a^lO^^GeV. 
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